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Abstract 

We study the dynamics of a quench-prepared domain wall state released into a system whose 
unitary time evolution is dictated by the Hamiltonian of the Heisenberg spin-1/2 gapped 
antiferromagnetic chain. Using exact wavefunctions and their overlaps with the domain wall 
state allows us to describe the release dynamics to high accuracy, up to the long-time limit, 
for finite as well as infinite systems. The results for the infinite system allow us to rigorously 
prove that the system in the gapped regime (A > 1) cannot thermalize in the strict sense. 

1 Introduction 

Statistical mechanics, be it for classical or quantum many-body systems, relies on the reasonable 
assumption that all microstates of a system in equilibrium are equally likely to be realized. In 
order to be correctly described by a statistical ensemble, a system must therefore be either coupled 
to an infinite thermal bath, or (if isolated) effectively display some form of ergodicity as per the 
eigenstate thermalization hypothesis d , 2 , thereby allowing it to forget its initial state by 'relaxing' 
towards a well-defined equilibrium state independent of the starting conditions. The existence of 
these relaxation processes is not usually questioned, but their precise mechanisms occur in a 
generally unspecified manner. For individual systems, classical mechanics distinguishes integrable 
and non-integrable cases, the former possibly having stable quasiperiodic orbits protected by the 
KAM theorem [3 , and therefore showing no true ergodicity. For non-integrable cases, ergodicity 
is taken for granted as a consequence of chaotic dynamics. 

In the quantum case, the question of equilibration has been the subject of many recent works, 
focusing in particular on the time evolution of a system after a quantum quench, at which a 
global parameter in the Hamiltonian is suddenly changed [1] - [53] . This is equivalent to releasing 
a prepared state and letting it evolve unitarily in time according to the Hamiltonian after the 
quench. Such situations, some of which are experimentally realizable, lead to many interesting 
theoretical questions. Under what circumstances is there a well-defined state at large times (i.e. 
what, if anything, remains of the initial conditions for given quench situations)? Can a general 
theory of relaxation based purely on the details of quantum dephasing be formulated? Are there 
different classes of systems {e.g. integrable versus non-integrable, finite versus infinite, gapped 
versus gapless) having different generic behaviour? What remains of the usefulness of commonly 
used quasiparticle bases to understand the dynamics after the quench? Clearly, much progress 
remains to be done in order to have a full understanding of all these issues. 

It is our purpose here to consider a well-defined situation, which can be handled to a sufficient 
degree of exactness to provide some partial but reliable insights into some of these questions. The 
approach we will consider is based on the use of exact wavefunctions for the Heisenberg magnet. 
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and therefore exploits the integrabihty of this system. On the one hand, this renders some of our 
results non-generic. On the other hand, we are able to provide hard facts for finite as well as 
infinite systems, gapped versus gapless, and exclude some possible scenarios. Since the theoretical 
understanding of the nonequilibrium dynamics of strongly-correlated systems is still in its infancy, 
such example cases will hopefully provide worthwhile reference points for later developments. 

The paper is organized as follows. In section 2, we first define our nonequilibrium problem, 
and present the tools we will use to study it. We then discuss in section 3 the set of eigenstates we 
use to obtain quantitative results. Section 4 considers the work probability distribution resulting 
from the quench for various system parameters. In section 5, we study relaxation dynamics by 
concentrating on the Loschmidt echo, which we compute. After a few words on the thermodynamic 
limit, we offer a discussion of our results and present our conclusions. All technical details for the 
computations are relegated to a series of appendices, which can safely be skipped by the reader 
only interested in the final results and conclusions. On the other hand, these appendices explain 
all calculations in sufficient detail to be reproduced by the specialist reader. 



2 Formulation 

We consider an isolated spin chain of N sites, each occupied by a local spin-1/2 degree of freedom. 
For definiteness, we put the spins on a ring and impose periodic boundary conditions. The thought 
experiment we perform consists in preparing the quantum state of the system at i = as 

|0)-|4_^t_^). (1) 

M N-M 

This state thus contains a magnetic domain wall between sites M and M + 1, and another (anti-) 
domain wall between sites N and 1 (in view of the periodic boundary conditions). This state 
can be prepared in different ways: we can view it as being created by an Ising model with an 
appropriate position-dependent field, or as resulting from a sudden polarizing pulse applied on a 
section of an initially fully polarized chain. For times t > 0, we let this state evolve unitarily in 
time under the antiferromagnetic XXZ Hamiltonian 



Hxxz — J 



N r 

E 



(2) 



Note that the Hamiltonian is rescaled by a factor of 1/A as compared to how it usually appears 
in the literature, in order the have a well-defined Ising limit (A — >■ oo). For A > 1 the spectrum 
of this theory is gapped, while for — 1 < A < 1 the system is in the quantum critical regime. 
We consider the case when the system is antiferromagnetic (i.e. J > 0). In order to simplify the 
expressions and without loss of generality we put J = 1 throughout the paper. 

The initial domain wall state ([I]) is not an eigenstate of the XXZ Hamiltonian away from the 
Ising limit A oo. On the other hand, for any fixed value of A, the exact eigenstates of this 
model form a basis in the Hilbert space on which we can at least in principle decompose the initial 
domain- wall state to arbitrary accuracy. Given this decomposition, the solution of the Schrodinger 
equation becomes straightforward, and we obtain the exact time-dependent wavefunction after the 
quench as the linear decomposition 



10(0) - J2 e^'^'-'Qnl*"). Qn = (^n|0), (3) 



f N \ . 

where the sum is over all the ( I eigenstates of Hxxz at fixed total magnetization, and the 

complex amplitudes Qn represent the overlaps (vector of the quench matrix) between the starting 
state and the exact eigenfunctions. Since we always work with normalized wavefunctions, the 



2 



coefEcients Qn should by definition satisfy tlie constraint 



(4) 



n 



Tliis constraint will constitute an important sum rule, allowing to quantify the accuracy of our 
results by assessing how faithfully the resulting wavefunction is reproduced. 

All the complexity of the problem is therefore hidden in two places. First, wavefunctions 
I'^n) and their energy En must be known. This is standardly handled by the Bethe Ansatz 
(see references [23] - [55] and references therein; we provide a summary of the necessary details 
in Appendices A and B). Second, the overlaps Qn must also be known. This is a problem of 
much greater complexity, which in the present situation finds its solution in the framework of 
the Algebraic Bethe Ansatz. The derivation of these overlaps is given in Appendices C and 
D. Note that similar overlaps (quench matrix entries) were also recently calculated using the 
Algebraic Bethe Ansatz for the case of the interaction-quenched Richardson model [501 HI] • Here, 
we also obtain these overlaps using integrability, but using a different method based on the explicit 
structure of the solutions to the Bethe equations. 

Despite being in possession of these two fundamental building blocks, one major difficulty 
remains. Since the Hilbert space is exponentially large in system size TV, the summation in (|3| 
is difficult to handle, and must in practice be truncated in order to reach sizes sufficiently large 
to allow the discrimination between finite- and infinite-size behaviour. We will show in the next 
section that this truncation is possible in the situation we consider. In particular, this puts us 
in position to study the long-time average of observables according to any prescription desired, a 
common one being 



and to make reliable observations on the relaxation/thermalization of the prepared state. 

The time evolution of this domain wall state has been studied in similar settings. For gapless 
spin chains this type of quench is studied in |29j . An exact analysis has been carried out for the 
XX-chain in [30]. The short time regime for the XXZ chain is studied using tDMRG in [3T], and 
an exact diagonalization was performed for the gapped XXZ chain with open boundary conditions 
[32] . We here offer a complement to these studies, consisting in results from integrability. 

3 Spectral analysis 

The space of eigenstates of the XXZ chain is spanned by Bethe wavefunctions, each described by a 
set of rapidity parameters {A} obtained as solution to the Bethe equations. Bethe wavefunctions 
are generically quite complicated objects, and can contain both unbound and bound magnons. We 
refer the reader to the standard literature [55J HT] [55] as well as to our appendices for the basic 
detail. In summary, solutions of the Bethe ansatz equations can be classified using strings (see 
Appendix [A] for a discussion). An n— string in a set of n complex rapidities sharing the same real 
part and invariant under complex conjugation. These strings can be interpreted as bound states 
of 'mass' n. In the Ising-limit it can be shown that an n-string corresponds to n adjacent down 
spins. From a perturbative analysis, we expect at large anisotropy A a well-defined hierarchy of 
states in terms of the importance of their overlap with the domain wall state. This hierarchy is 
presented in Table [l] The most important states are those consisting of a single M— string, which 
occupy a dispersionless line in the thermodynamic limit. They are followed by successively more 
complicated partitionings of the M rapidities into more and more individual bound or unbound 
states, each partitioning representing a whole continuum of excitations. The number of particles 
in the state is defined here to be the number of elements in the partitioning. From a perturbative 
analysis it can be easily shown that for sufficiently large M the system is insensitive to small 
changes of M {i.e., of how distant the domain and anti-domain walls are from each other). This 
allows us to perform the calculation with an M slightly below N/2 without loss of generality. 




(5) 
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Order (1/A) 


strings 


^ particles 


Energy (Ising limit) 





{M} 


1 


-1 


1 


{1,M-1} 


2 


-2 


2 
2 


{l2,M-2} 
{2,M-2} 


3 
2 


-3 
-2 


3 
3 
3 


{l3,M-3} 
{l,2,Af-3} 
{3, M - 3} 


4 
3 
2 


-4 
-3 
-2 











Table 1: The hierarchy of states, in order of importance of contributions to the normalization sum 
rule Q. Each partitioning of M rapidities into strings leads to an independent excitation class, 
only the few shown here being of relevance to the current setup. 



The reason for this choice is that for these values the string hypothesis is better satisfied for 
the majority of string states. In the left panel of figure [T] we give the band structure of the 
dominant classes of excitations as a function of A. The number of particles corresponds to the 
number of strings a state is made of. The right panel of the same figure shows the normalization 
sum rule saturation coming from these excitations. For sufficiently large values of anisotropy, the 
saturation is essentially complete using these states only. Approaching the gapless regime however 
requires increasingly large numbers of classes of excitations, since all overlaps scale to zero as 
can be expected from Anderson's orthogonality catastrophe scenario |33j . A further issue with 
the gapless limit A 1 is that in its vicinity, solutions to Bethe equations involving complex 
rapidities become more and more difficult to find, and often require considering deviated strings 
in detail '34]. We will not burden ourselves with these issues here, and only consider systems 
sufficiently deep in the gapped regime where such deviations are negligible. 




A 



Figure 1: Left: sketch of the energy values covered by the highest three energy bands of the XXZ 
model as a function of A compared with the energy of the initial state Eq. Indicated is the number 
of particles the states in different bands are made of. Right: normalization sum rule saturation 
coming from summing over all states of each basic excitation family, as a function of A. As the 
gapless regime is approached, all overlaps scale to zero. 



4 



4 Work probability distribution 



One of the most straightforward measurable quantities which can be obtained from the knowledge 
of the overlaps Qn is the work probability distribution 35J defined as: 

n 

The work expectation value can easily be computed analytically: {W) — {<j>\Hxxz ~ -ffol^) — 
and {W'^) — 232- • (W) can thus be considered intensive here, which for a generic quantum quench 
is not the case. In figure [2] we plot P{W) for two different values of A. The calculations are done 
for a finite system size {N = 250, M = 100), therefore the results are binned in energy (using bins 
larger than the interlevel spacing) in order to generate smooth curves. The contribution of the 




Figure 2: Work probability distribution for N — 250 and M — 100 and anisotropies A = 1.5 (left) 
and 3 (right). Indicated are the contributions from the important string states. The vertical scale 
of the full peak of the M— string is cut in the plot for clarity. 



M— strings is a narrow peak with a vanishing bandwidth in the thermodynamic limit, see (18). 
The 1, M — 1 states have two peaks at the edge of the energy band. The peak at the right of the 
continuum close to {W) is the result of individual peaks associated to a large value of overlap. The 
peak at the left of the continuum is the result of smaller overlaps accompanied by an increasing 
density of states. Other states made up of two strings lie within the energy band of the 1, M — 1- 
strings and display a similar structure. Since the energies are intensive in this energy domain it is 
expected that the results for P{W) we obtain closely mimics the one in the thermodynamic limit. 
This is motivated by considering the difference in P{W) for two different system sizes, see figure 
[3] showing the fact that finite-size effects rapidly disappear. 



5 Relaxation dynamics: Loschmidt echo 

We now turn to the question of whether the system effectively displays some form of relaxation. 
Of course here, since we consider a single realization and not an ensemble in the presence of a 
bath, relaxation can only come from the relative dephasing of the various terms in the right-hand 
side of (|3|. 

It is not a priori easy to guess the outcome of the time evolution, since the states occupy 
coherent modes followed by overlapping continua, and each state contributes according to its 
relative overlap. Moreover, this outcome can essentially depend on which observable is considered, 
via the form factor values of the operator concerned, each operator favouring specific sets of state 
combinations. As a concrete example, we will focus on studying the Loschmidt echo [35l|36j which 
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A=1.5 A=1.5 




-1.5 -1.0 -0.5 0.0 -1.5 -1.0 -0.5 0.0 

W W 



Figure 3: For A = 1.5 we plot the difference of P{W) obtained from different system sizes. Left: 
the difference between N = 200 and N — 150. Right: the difference between N = 250 and 
N — 200. The differences fall off with increasing system size, and for the sizes presented are 
smaller by about two orders of magnitude than the value of P{W) obtained (see figure [2]). 



is the overlap between the inital state and the time evolved state, 

m,n 

2 

(7) 

This will allow us to quantify if the system dephases, if so how quickly and to what point, and 
whether initial state (partial) revival can take place. Plots of the Loschmidt echo for short and long 
times are presented in figure[4] We can separate different regimes in the time dependence. First and 
foremost, the small time, transient regime plotted in the left panel can be easily understood using 
simple perturbation theory, which predicts an initial decay quadratic in time with an anisotropy- 
dependent coefficient: 

C{t) ^ 1 + {{Hf - {H^))t^ + 0{t^) 

= ^-^t' + 0{t% (8) 

More complex behaviour becomes evident when considering longer times, as in the right panel 
of figure |4] One point worth emphasizing is that our method, unlike e.g. real-time numerics, is 
(in view of the exact wavefunctions and energies we use) directly applicable irrespective of the 
time t considered, and thus valid also at large times. There however, finite-size effects show up, 
originating from the discrete nature of the energy levels. Numerical inaccuracies also prevent 
us from giving numerical values at arbitrarily large times. Returning to the figure, a number of 
interesting features are displayed. First, interferences between the various terms lead to oscillations 
around the long-time average 

£ = J2S{Em - S„)|(0|*„)n(0|*,„)|2. (9) 

These oscillations are accompanied by a decaying envelope which can be interpreted as effective 
relaxation. This behaviour carries on until a new regime is entered, where finite-size effects take 
over and lead to accidental partial revivals where the Loschmidt echo shows rapid, larger-scale 
variations. More generally, the echo looks chaotic after the onset of these revivals, and the numer- 
ical calculations lose their meaningfulness. The revival time is observed to grow with system size, 
first approximately linearly with system size, but more generally in an unspecifiable and nontrivial 
fashion. 



J P{W)e'^'dW 
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Figure 4: Loschmidt echo as a function of time for the domain waU quench for A = 1.5 and 
N = 200, M = 80. In the left panel, the short time behaviour is shown. The dashed line is the 
result from perturbation theory. At large times, plotted in the right panel, one can distinguish 
effective relaxation towards the long-time average, but also finite-size induced small-scale revivals 
after a certain point (see main text). Note the change of scale in the horizontal axis in going from 
one figure to the other. 



6 Thermodynamic limit 

The overlaps of the M— string states with the initial domain wall state can be computed when 
N oo both for the case of M finite and for M — )■ oo. When M oo the M-strings become 
exactly coherent. We can compute a lower bound for the long time average of the Loschmidt echo 
C. In appendix D we derive: 



N 

hm ^|(</'|OP = n(l-^"'"'')'' ^>1' (10) 



in which we have used the definition cosh ry = A, and where the summation on the left-hand side 
is taken over all Af-string solutions. By only taking these contributions into account we obtain a 
lower bound for ([9]): 

oo 
ri=l 

This result can be interpreted as the overlap between the initial state and the asymptotic state. 
The non-zero overlap indicates that the asymptotic state keeps most of the spatial anisotropy from 
the initial state. We can conclude that for this type of quench the system will not thermalize in 
the thermodynamic limit. 

In table [2j we compare this result for C with the one for a finite system with N = 200. The 
finite time average is over intervals before the finite size effects show up. The finite size result 
-^200 exceeds the thermodynamic limit value Coo by a small value, which can be attributed to the 
transient regime, which can be neglected in the infinite time average. This clearly suggests that 



the lower bound ( 11 ) is the exact value for the infinite time average in the thermodynamic limit. 



7 Discussion and Conclusions 

In this work we studied the relaxation dynamics of the gapped XXZ chain after a quench from 
an initial domain wall state. The techniques that were used are based on the Algebraic Bethe 
Ansatz for a finite chain. Some of the results were extended to the thermodynamic limit. The key 
ingredient for studying this quench was deriving a numerically efficient expression for the overlap 
between the initial state and eigenstates of the XXZ Hamiltonian. Overlaps for different Ising 
states as initial state could be computed following the same logic, however in the general case of 
a finite density of domain walls this will become an intractable combinatorial problem. 
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A 


c 


^200 


1.5 


0.481212 


0.481402 


2 


0.725941 


0.726285 


3 


0.884184 


0.884241 


4 


0.936021 


0.936077 



Table 2: The lower bound for C in the thermodynamic hmit compared with long time average for 
a finite size system {N 200 and AI = 80). 



Since the energies of the relevant part of the spectrum are intensive a good comparison with the 
thermodynamic limit is possible. One important result is that the bandwidth of the M— strings 
is nonzero for finite systems allowing the long time average of the Loschmidt echo going to zero, 
something that does not happen in the thermodynamic limit. Hence we conclude that the limits 
limr_>.oo and limAr_>.oo do not commute, which was also pointed out in |36j . 

In the thermodynamic limit we derived a lower bound for the infinite time average of the 
Loschmidt echo and argue that this bound is the exact value. This led us to the most important 
conclusion of this paper, namely that the system for this quench will not thermalize. The reason 
why the system does not thermalize can be explained from its spectrum, and does not directly 
rely on the fact that the system is integrable. Imagine for instance a gapped system with lowest 
(or highest) band having zero bandwidth in the thermodynamic limit, such as a system with 
macroscopically-degenerate ground states. Consider a quench such that (W) lies right between 
the lowest (or highest) energy band and the next one. A lower bound for the overlap K^jvPo)!^ 
of the lowest band can then be obtained by simple reasoning. Let Eq be the energy of the lowest 
energy band of zero bandwidth. We denote the band gap by Ei — Eq. The expectation value of 
the initial Hamiltonian is (Hq) = Ef) + AE. We then choose < AE < Ei — Eq. In this case we 
can obtain a lower bound for the total overlap of the Eq states. 

£;o + Ai; = ^|((^|^'„)p£^n 

n 

> |(./)|^'o)|'£;o + (i- K'/'l^'o)!')^^!, (12) 

from which follows the lower bound 

|(<^|*o)P>l-^;^ ^. (13) 

In other words, such a quench must occupy the lowest (resp. highest) dispersionless band with 
0(1) amplitude. Since these states are dispersionless and separated from other states by a band 
gap, they cannot dephase, and therefore cannot effectively relax. The large-time asymptotic state 
will therefore always maintain memory of the ground state occupation amplitude of this band, 
which is a quench-dependent quantity. Note that such situations trivially cannot show eigenstate 
thermalization according to the arguments of in which such circumstances were argued to be 
non-generic and therefore not treated. 

The condition used for deriving this lower bound is quite strong. In figure [T] one can see that 
for 1 < A < 1.5 that {W) does not lie in the gap anymore, but the overlaps of the M— string 
are still substantially big. One can also question whether the condition of an isolated peak is 
essential. From the work probability distribution it is clear that the contributions of states made 
up of different strings can be distinguished. Simply looking at energy levels is therefore insufficient 
to assess the presence or absence of relaxation, and one must also take into account the values 
of the wavef unction overlaps. If we consider for example n strings with a diverging length in the 
thermodynamic limit, they all will become coherent. Although they are embedded in a continuum, 
we can reasonably expect a non-zero height peak of zero width in P[W), which again will prevent 
thermalization. 
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In the case of evolution under a gapless Haniiltonian, our method loses its efficiency, since all 
overlaps scale to zero and the number of states to take into account grows accordingly, preventing 
an efficient truncation of the sum in |3] The absence of a gap thus opens the door to dephasing 
involving arbitrarily complex excitation continua, in correspondence with Anderson's orthogonality 
catastrophe principle. Therefore it is expected that £ = in this case, and the possibility of 
thermalization remains, although it cannot be quantified here. The presence of a gap is thus 
determinantal to the long-time asymptotics of dynamics after the quench release of the domain 
wall state, but is not very sensitive to system size for large enough systems. 

To summarize we considered a quench starting from a domain wall state under the time evolu- 
tion of the gapped XXZ chain. We observed that the long time behavior cannot be described by 
a statistical ensemble. This can be understood by the fact that the spectrum is gapped, and does 
not directly rely on the fact that the system is integrable. Furthermore we were able to assess finite 
size effects by considering both results for the finite and infinite systems. General conclusions for 
the long time behavior of the XXZ model after quench cannot however be drawn from the results 
presented here. For instance, the two quasi-degenerate ground states of the XXZ chain we are 
dealing with are separated from all other states by the gap, so one might expect a similar effect if 
one would perform a quench starting from the Neel state. However, the low energies are extensive 
in contrast with the high energy excitations which are intensive. The effect of the gap is much 
less pronounced and one might at least expect that the staggered magnetization will vanish for 
all A > 1, in correspondence with the results of [TO]. In future work, we will consider this kind of 
initial state as well as the time evolution of more elaborate observables. 



A Bethe ansatz equations 

In the Algebraic Bethe Ansatz the eigenstates are completely characterized by a set of rapidities 
{Xj},j = 1 . . . M, solutions to the Bethe equations. These solutions can be either real or complex. 
The complex solutions typically come in groups, invariant under complex conjugation, which 
represent bound states. At a small density of down spins or large anisotropy A, most solutions 
take the form of so-called strings, in which a number of rapidities share a real centre while the 
imaginary parts are equally spaced according to the so-called string hypothesis. An n— string takes 
the form: 



XZ + if]{n + 1 ~ 2j)/2 + i5j 3 



(14) 



where r] = acosh(A), and 5j is a deviation. When the string hypothesis holds (i.e. when all 
deviations 6j are sufficiently small) the Bethe equations can be rewritten as equations for the real 



parts of the rapidities only, yielding the Bethe-Gaudin-Takahashi equations 
version of the Bethe-Gaudin-Takahashi equations for A > 1 is given by 



.(a:^) 



27r 



1 

N 



E 0^ 

(m,/J)#(«,Q) 



- A^) 



The logarithmic 



(15) 



with the dispersion and scattering kernels 

^ri(A) = 2 arctan 



tan A 



(^tanh 71,77/2 



■2ti 



0mn(A) 



(A) 



71 7^ ni 
n = ni. 



(16) 



(A) + 2e\ (A) + . . . + 26'„+„i_2(A) 

^2^2 (A) -f 204(A) + ... + 202„-2(A) + er,+m{\) 

Solutions of the Bethe-Gaudin-Takahashi equations can be classified by the set of quantum num- 
bers which are integers (half-odd integers) if TV — Mn is odd (even). The partitioning into 



strings satisfies the constraint 'Ylin=i = M. The energy of an n— string A^ is 

sinh nrj 



e„(AS) = -tanh(77) 



cosh 7177 ~ cos 2A2 ' 



(17) 
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The total energy and momentum of an eigenstate containing strings are given by 

P = ^^Af„-^^/^ mod 271. (19) 

n a.n 

For finite rj and large n the energy e„(A) becomes dispersionless. An n— string can be interpreted 
as a magnon bound state with a mass proportional to n. 



Band structure 



In the thermodynamic limit {N oo) the band structure is easily obtained from (18). For a 
given string configuration {Af„} one can define: i?™"^({M„}) by putting all AJJ to tt/2. Similarly 
one can define E™-^"^ {{Mn}) by putting all AJJ to 0. The M— string is a special case since it has a 
vanishing bandwidth in the limit M ^ oo. In the Ising limit (77 — >■ oo) the energy of a state made 
of m different strings is —m. This corresponds in the Ising language to a state made of m distinct 
blocks of down spins and clarifies the bound states interpretation in this limit. 



B Classification of string solutions 

In order to evaluate the spectrum numerically a correspondence between quantum numbers and 
rapidities is needed. For the isotropic chain (A = 1) every set of quantum numbers corresponds to 
a unique solution of the Bethe-Gaudin-Takahashi equations. A complete classification of solutions 
in terms of quantum numbers was given in [28 . For the case of the XXZ chain with A > 1 
this classification in not known in general. For solutions with only real solutions this was done 
in [10]. Here we generalize the result to string solutions. Note that a genuine proof that the 
Bethe-Gaudin-Takahashi equations produce the right number of eigenstates for A > 1 remains to 
be found. 

Bandwidth: 

In |3D] it was argued that Xj — Xk if Ij < Ik- Equal rapidities do not yield proper Bethe wave- 
functions, and we therefore only need to consider sets of distinct quantum numbers. In the 
parameterization chosen in the Bethe-Gaudin-Takahashi equations all rapidities are restricted to 
an interval of width tt, i.e. for all strings of given length n we have: A^ — A" < tt. From this 
requirement and assuming that the Bethe-Gaudin-Takahashi equations are monotonic in A^ it 
follows that: I^j — /" < iV — J2m ^nniMm = Wn with = 2min(ri, m) — 6n,m- We call Wn the 
bandwidth of the quantum numbers {If}- 

Quasi-periodicity of the Bethe-Gaudin-Takahashi equations: B„ 

Since A + tt = tt as far as the wavefunctions are concerned, there is no one-to-one mapping of 
quantum numbers and wavefunctions. Namely, the shift A" — ^ A" -I- tt leaves the wavefunction 
invariant although the quantum numbers in the Bethe-Gaudin-Takahashi equations change. For 
every string length n we can define a transformation: S"" : {(A°,/°)} — ?► {(A°,/°)} with: 

j = 1 . . . M„ - 1 
m 7^ n, j = 1 . . .M„. (20) 

From this transformation we see that /">/"-!- 2n. For every string length n there are Bn = 2n 
bands of quantum numbers which cannot be transformed into each other by 5", whereas in the 



(Ki^JmJ - (Ar+^,/r-H(I^„ + (2n-l))) 

(a;',/;) = (a^i,/;Vi+(2^-i)) 

(AJ*, IT) = i^T ^ "I" tnm) 
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XXX case there is only one such band for every string length. If we restrict to only one string 
length it follows that the number of choices for the quantum numbers {^"} in a single band is: 

. In case of Bn bands the number of solutions is: 



Mr 



Degeneracies from different string sectors 

To find the minimal set of quantum numbers one should first determine using 5*" the bandwidth 
Wn and the number of bands i?„ for every string length n. However, the transformation S*" not 
only affects the quantum numbers of length n but also all the others. To take this effect into 
account one should also consider transformations of the form S''"(S'")~^, leading to additional 
state exclusions. In general these are very difficult to determine. The explicit classification for 
two string types is given below. 



Two types of strings 

We consider m n— strings and m n— strings such that: mn + fhn = M. The bandwidths of the 
two sets of quantum numbers are: 

Wn ^ N - {2n - l)m - 2nm (22) 
Wn ^ N- 2nm - {2n - l)m. (23) 

The number of bands are _B„ = 2n and Bn = 2n. To take care of the degeneracies coming from 
different string sectors we consider the transformation: 



n^^yi ri-l...I{^ + Wn-l \ 

l/« -W-ft-(2n-l) + 2n.../,|_i-(2n-l) + 2n; • ^ > 



This turns out to be the only extra transformation one needs to consider. Define: I^ax — 
{Wn-l) + {Bn-l). From the conditions: + VK„ - 1 < I^^^x and -Wn-{2n-l) + 2n> /, 
follows that we have to exclude: e {F^^^^ . . . I^^^ + 2n - 1} and S {/^^^ - 2n + 1 . . . I^^J 
The total number of unique solutions for this case is 



min ' 
n 

min 



{Bn - l)m + Wn fWn\ {Bn - + Wn (Wn\ \2 f Wn ' l\ f Wn - 1 

( 1 _ (2n) 



Wn \ m J Wn \ m J Km—llKm — 1 

- 2n{m + m)7V (Wn\ f Wn 



W„Wn \ m J \ m 



(25) 



C The overlap matrix 

The Algebraic Bethe Ansatz deals with the problem of diagonalizing simultaneously the transfer 
matrix T(A) for all values of A, from which all conserved quantities can be obtained [2^. All 
operators can be expressed in terms of four non-local Hilbert space operators A{X),B{X),C{X) 
and D{X). The Algebraic Bethe Ansatz requires a pseudo- vacuum |0) which in case of the XXZ 
chain is the state with all spins up, such that: 

A(A)|0)=a(A)|0) 
D(A)|0)=d(A)|0) 
B(A)|0)^0 

C(A)|0)=0 (26) 
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The transfer matrix is constructed as: T(A) = {A + D){X). The vacuum eigenvahies of the A{X) 
and -D(A) operators are: 



A + i A = 1 

</5(A + i??) = <( sinh(A + i?7) |A = acos(77)| < 1 (28) 
^sin(A + 11]) |A = acosh(ry)| > 1. 

The set are inhomogeneity parameters, which should be set to S,j = if]/ 2 in order to obtain the 
XXZ chain. States, respectively dual states can be constructed as: 

M 



M 

= {Q\X{C{\,). (29) 



In order to represent an eigenstate the rapidities \j should satisfy the Bethe equations. These 
states are not automatically normalized, we therefore write: \^) = IV') / V (V'lV')- The norm of a 
Bethe state in case of strings is given by [261 EZl l38l l39] : 



{m = ^{^vr n ^\+;^^ det<i>(-)({A,}) + 0(J) (30) 

where <i>*^'') is called the reduced Gaudin matrix: 

= s.kS^P i^:^^^^^^^ - E - a!,) I +(i-<5,fc<5„,)-^e,fe(A^„-A^). 



(31) 



An extremely useful formula is Slavnov's expression for the scalar product between an eigenstate 
(■^({A})! (represented by a set of rapidities {A} that satisfy the Bethe equations) and an arbitrary 
state \tp{{lJ,})) (without restrictions on {/i}) JH]: 



with 

ip{iT]) 



/' M M \ 

W fih - 6 + iv) - d{]ib) Y\ 'Pi^i - f^b - iv) I • (33) 
^/=i 1=1 / 



C.l General overlap 

The initial state can be written in the Algebraic Bethe Ansatz language using the inverse mapping 

j-l N 

Si ^Yl{A + D)iCk)B{(,) YliA + om- (34) 

fc = l 1=3 + 1 
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Making use of n^i(^ + = 1 and {A + D){^j)\0) = |0) we can write the state as: 

M 

l'^) = n^(^^-)|0) (35) 

which is a normahzed state. The overlap between an eigenstate for finite A and the initial state 
\(f)) can be computed using Slavnov's theorem for the scalar product. In the case when {/ij} — 
the matrix H simplifies: 

Hat = -T , ,'^^,7^ , ^ . , n ^(^' - + (36) 

For the XXZ spin chain all inhomogeneity parameters should be set to = ir]/2 V j. The limit 
should be taken carefully, there are M{M — 1) poles coming from: Iljxfe '/'(O ~ ?fe)- Since the 
columns of the matrix H are also M(Af — l)-fold degenerate we can apply I'Hopital's rule. First we 
notice that we can easily extract factors Jli=i vi^i ~^b + i'n) from the determinant, since applying 
I'Hopital's rule to those terms does not remove any degeneracies. 

^ f(M . (38) 

Note that this determinant of H is the same as that in the partition function of the 6- vertex model 
with domain wall boundary conditions [13]. However in this case both the sets {Xj} and {^j} play 
the role of boundary conditions and satisfy no Bethe equations. To be able to apply I'Hopital we 
are interested in 9? Hab- First we rewrite Hab- 



</5(Aa - Cb) "^(Aq - 6 + "i-V) 



Next, consider 



<^(Aa - 6) / V <^(Aa - 6) / V '/'(Aa - 6 

Using this result it is easy to derive for n even: 

-U(A.-eb + «;)j I' ^''^ 

and for n odd we get a similar result. In general the coefficients c" are complicated expressions. 
However it is straightforward to derive that = ^O'' even. Since in this case there is no 
mixing between the left and right terms in (40). It will turn out that theses coefficients are the 
only ones we need, since all terms with coefficient c" for j < n/2 can be removed using row manip- 
ulations while leaving the determinant invariant. Differentiating the denominator: Ilixfc '/'(Cj 



k 



is much simpler. First we put — irj/2. If we now consider ^2 we get a single zero in the limit 
^2 = *^/2. We can continue this logic and see that gives rises to j — 1 zero's, hence we get a 
factor — We also notice that the factors n\ from the numerator cancel against the factors 
from the denominator so we are left with: 



detiJ 



m) =ny(A/+V2)"-^ ^ "7f , , (42) 
f^i nj>fc<y5(Aj - Afe) 



I/.b=(4^^V-f4^^V- (43) 



^{Xa-lll/i)) \^{\a+iTl/'i) ) 
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The matrix H has a similar structure to a Vandermonde determinant. This makes the matrix 
iU-conditioned and is therefore not suitable for a numerical evaluation. However, for a given string 
structure the determinant can be written in terms of sums of Vandermonde determinants. In the 
following sections this is done explicitly for the most important cases. The normalized overlap is: 

{m Y{^^,^{\l+^^/2)^AeiH ^^^^ 



VW^ y^_lM(Af-l)/2^(i^)A/ Y{._^^^ ip{Xj - Afe + if]) det $ 



These expressions hold for all three paramcterizations in ( 28 1 



C.2 Overlap in case of an M— string 

In the case where the solution is an M— string: X^j = Aq + ir]{M + 1 — 2j)/2 the determinant of 

AVan ' 

all Xj distinct. Its determinant is given by: 



H is exactly a Vandermonde determinant. A Vandermonde matrix is defined as: Vjk = ^ with 



detV = Ylix,-Xk). (45) 

3<k 



Introducing the notation 

aj = ; j = i...m + i, (46) 



^(A,+z7?/2)' 



the determinant of H can be written as an Af + 1 x Af + 1 Vandermonde matrix determinant: 
det H = det 




_M / \ 1 „ „M 

yClM+i — aM ■■■ 0,M+1~^M/ \^ O.M+1 ■■■ O^M+l/ 

*ff' ( ifi'jX^' + iviM + 2~2j)/2) _ y'(A^^ + »7?(M + 2-2fc)/2) \ 
^\ \ Vi^^ + MM + 2 - 2j)/2) ^(Af + iriiM + 2 - 2A:)/2) J 



The reduced Gaudin determinant for an M-string takes the simple form: 



(47) 



,M _ M <p{ir]M) 

ifix^i - iM'n/2)^{x^^ + iM'n/2) ■ 



det = ..ZiZ^Im , ...^/ox - (48) 



This results in the normalized overlap for an M-string: 



{^\^) _ nf=iV(*"'7) /^(A„-ir7M/2)^(A„ + i?7M/2) 



^(A„ - ir7M/2)Af V _iM(M-i)/2^ 
The absolute value squared of the overlap is 



(49) 



(50) 



C.3 Overlap for two (m, m)-strings 

We consider states constructed from two strings with lengths m and fh such that m + m = M, 
and m <fh: A™ = A™ + iT]{m + 1 - 2j)/2 and Xf = A" + irjiffi + 1 - 2fc)/2. Some convenient 
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notations are: 



_ (^'(A;" + ir7/2) 



Now the overlap matrix takes the form: 



/ a2 - ai ... af - af \ 



det 



62 - hx 



''2 



"m+1 "m / 



7 = 1 ... TO + 1 
j = 1 . . . TO + 1. 

/I 1 ai 



det 



1 1 flm+l 

1 61 
\0 1 



M 



(51) 
(52) 



(53) 



Expanding this determinant in the first cohimn results in a sum over m + 1 Vandermonde deter- 
minants. After factoring out common factors we obtain: 



m+l / m+1 m+1 



m+1 



m+l 



det 5 = ^ [] n ("'^ - n («p -"9)1 n ~ ^-)- ^^i) 
i=i I 1=1 



P < Q 
P, Q ^ j 



C.4 Overlap for a general string solution 

Consider a state constructed of n strings with lengths: li < I2 ■ • ■ ^ ^n- Introducing the notation 



^^(Ai°^.+z77/2) 



j — 1 . . Aq, + \ a — 1 . . .n, 



the overlaps can be written as a sum over Vandermonde determinants: 



(55) 



deti/=E--- E n n nK-«f)n n (^e) 

jl = l j,i-l = l \Q!<;3 j^jc fc/j^) 1=1 j<k j.k^j^ 



D Thermodynamic limit 

The Bethe-Gaudin-Takahashi equations for the M— strings reduces to the simple equation: 

n _ j-M 



(57) 



For even M and A > 1 the A'' possible quantum numbers are {—{N — l)/2,—{N — l)/2 + 



1, . . . , (A*" — l)/2} which follows from (20 1. In the thermodynamic limit we can define a density 



function p{X) for an M-string solution centered around A: 
n(\'') ft - ^ 



-iif{i'qM) 



27r d\^^ 



27r ^(A^'^ - i'qM/2)ip{\M + iriM/i) 



(58) 



such that p('^)'^'^ — This function p(A) should be interpreted as the number of solutions 

there are with a rapidity in the domain A, A + dX. The total contribution of all A^ M— string 
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solutions is: 



p{X)dX 



2Fi(-,M,l,^'(zr?M/2)-2) (59) 



where 2-^1(0,6,0,2;) is the hypergeometric function: 



2fi(a,6,c,z) = ^ (rc)„ = a;(a; + l)(x + 2)...(x + n-l). (60) 

In the limit when both M and A'' are sent to infinity a simpler expression can be obtained: 



M,N^oo^^ i^j^l^j^) M^oo J_^/2 ((/?(ir?M)/2)^-i Trip'{ir]M) \ip'{ir]M) 
= lim y +0^ ^ 

M^oo IJL (^'(i,yM)/2 



n (61) 

Tl=l 



Thermodynamic limit: XXX 

The isotropic limit A 1 should be taken with care, since some M-strings solutions correspond to 
difl[erent string configurations. For the XXX case we have the following restriction for the quantum 
number of an M-string: |/^| < {N — M)/2. Prom the Bethe-Gaudin-Takahashi equation: 



7M 



(A^) = 27r7^/7V (62) 



we see that in the limit N ^ 00 and finite M allowed values for A go from minus to plus infinity. 
In this limit we can derive the weight function: 

^-2^dA^^^(^ 2;^(AM-iM/2)(A^+iM/2)' ^^^^ 
so the contribution of the M-string for finite M becomes: 

''^'Uii't'm' r.x (M-i)p 



\w'm' ^ r (M-i)p 



WW) 

_ (2M- 2)! 
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